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SECTION  1 


INTRODUCTION 

Elasto-plastic  models  with  non-associated  flow  rules  are  commonly  used  today  in 
the  defense  community  to  describe  the  responses  of  soils,  and  other  complex  geologic 
media,  to  explosive  sources.  Originally,  non-associated  flow  laws  were  introduced  into 
such  models  for  the  purpose  of  bringing  the  predicted  volume  changes  more  in  line 
with  observed  changes.  However,  in  the  past  several  years,  it  has  become 
increasingly  apparent  that  there  are  more  compelling  reasons  for  using  non-associated 
flow  laws  in  conjunction  with  soils.  First,  there  is  considerable  experimental  data  now 
which  show  that  plastic  flow  of  granular  materials  near  failure  is  not  associated 
(Drescher  and  de  Josselin  de  Jong.  1982:  Oda  and  Konishi,  1974:  Drescher,  1976: 
Lade,  Nelson,  and  Ito.  1987)  and,  secondly,  it  has  been  theoretically  established  that 
frictional  materials,  in  general,  exhibit  non-associated  flow  laws  (Mandel,  1947; 
Spencer.  1964:  Rice,  1976).  Inasmuch  as  soils  are  frictional  materials,  in  which 
deformation  occurs  through  the  frictional  sliding  of  granular  particles,  there  appears  to 
be  both  a  theoretical  and  physical  basis  for  adopting  non-associated  flow  laws  in 
modeling  soil  behavior. 

In  granular  materials  the  phenomena  which  invalidate  the  adoption  of  associated 
flow  rules  are  termed  "non-normality"  and  "non-coaxiality".  Non-normality  implies  that 
plastic  strain  rate  is  not  normal  to  the  current  yield  surface.  Non-coaxiality  means 
that  the  principal  axes  of  stress  and  plastic  strain  rate  do  not  coincide,  which  usually 
leads  to  the  non-normality.  An  extensive  review  of  microscopic  and  macroscopic 
plasticity  theories  for  geomaterials,  which  exhibit  non-normality  and  non-coaxiality  was 
compiled  by  Nemat-Nasser  (1983).  In  phenomenological  plasticity  theories  the 
noncoaxiality  has  been  introduced  either  by  the  yield-vertex  model  (Rudnicki  and  Rice, 
1975)  or  by  a  smooth  yield  surface  model  with  kinematic  hardening.  In  this  report, 
only  plastic  strain  rates  expressible  by  the  flow  potential  are  considered  together  with 
isotropic  hardening  of  smooth  yield  surfaces.  Therefore,  the  component  of  plastic 
strain  rate  which  represents  the  effect  of  non-coaxiality  without  flow  potential  is  not 
treated. 

Recently,  non-associated  flow  laws  have  been  the  subject  of  considerable 
discussion  and  debate  in  the  ground  motion  community.  Specifically,  non-associated 
flow  laws  have  been  attacked  on  the  basis  that  they  violate  ll’iushin’s  Postulate  and 


1 


lead  to  non-uniqueness  of  initial  value  problems.  Since  ll’iushin's  Postulate  is  not 
grounded  in  the  laws  of  thermodynamics,  as  we  shall  show  in  the  sequel,  its  violation 
is  clearly  not  a  critical  issue.  The  question  of  uniqueness,  however,  is  of  concern. 

The  purpose  of  this  report  is  to  explore  the  general  restrictions  that  are  placed 
on  non-associated  flow  laws  by  the  second  law  of  thermodynamics  and  by  the 
requirement  that  unique  solutions  exist  to  initial  value  problems.  The  basic  question 
that  we  seek  to  address  is  as  follows:  Are  there  any  restrictions  imposed  by  either 
thermodynamics  or  uniqueness  which  make  non-associated  flow  laws  generally 
unacceptable  for  use  in  conjunction  with  explosive  ground  motion  studies? 

In  Section  2.  the  basic  relations  from  the  Onsagerist  approach  to  the 
thermodynamics  of  irreversible  media  are  developed.  The  postulates  of  Drucker  and 
ll'iushin  are  discussed  in  Section  3,  where  it  is  shown  that  they  are  unrelated  to  the 
second  law.  In  Section  4,  the  basic  features  of  non-associated  elasto-plastic  models 
are  explored,  including  continuity,  path  dependence  and  a  peculiar  feature  in  which  the 
plastic  tangent  modulus  can  exceed  the  corresponding  elastic  modulus.  The  question 
of  uniqueness  of  initial  value  problems  for  non-associated  elasto-plastic  solids  is 
examined  in  Section  5.  where  it  is  shown  that  the  usual  approach  for  proving 
uniqueness  is  unsatisfactory  when  the  plastic  flow  law  is  non-associated.  Finally,  in 
Section  6.  two  classes  of  non-associated  elasto-plastic  models  that  are  presently  being 
used  by  the  defense  community  are  considered,  and  the  implications  of  the  general 
restrictions  imposed  by  thermodynamics  are  discussed. 
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SECTION  2 


THERMODYNAMICS  OF  IRREVERSIBLE  MEDIA 

In  this  section,  the  restrictions  imposed  by  thermodynamics  on  elasto-plastic 
constitutive  theories  having  non-associated  flow  laws  are  explored.  The  two  main 
approaches  to  the  thermodynamics  of  irreversible  media  are  briefly  reviewed  and 
difficulties  are  noted.  The  approach  currently  accepted  by  the  majority  of  researchers 
in  the  field,  namely,  the  Onsagerist  approach,  is  adopted  in  the  present  study  and 
forms  the  basis  for  the  developments  which  follow. 

2.1  APPROACHES. 

There  are  essentially  two  major  schools  of  thought  that  have  emerged  in  recent 
years  regarding  the  thermodynamics  of  continuous  media,  namely,  the  Onsagerist 
irreversible  thermodynamics  of  internal  variables  (Onsager.  1931;  Eckart.  1948;  Meixner, 
1953;  Biot.  1954;  Zeigler,  1958;  Schapery.  1964:  Valanis,  1968,  1971c;  Rice.  1971; 
Nemat-Nasser.  1975a, b)  and  the  rational  thermodynamics  (Truesdell  and  Toupin.  1960; 
Coleman,  1964;  Truesdell  and  Noll,  1965).  The  essential  differences  between  these 
two  approaches,  which  have  been  discussed  by  Kestin  and  Rice  (1970),  Valanis 
(1971c)  and  Nemat-Nasser  (1975b),  lie  in  the  concepts  of  entropy  and  temperature. 
As  a  result,  they  lead  to  different  inequalities  implied  by  the  second  law  of 
thermodynamics  (Valanis,  1971c). 

In  the  rational  thermodynamics  approach,  temperature  and  entropy  are  assumed 
to  exist  strictly  as  mathematical  entities,  without  physical  interpretation,  and  a 
Clausius-Duhem  inequality  is  postulated  by  introducing  suitable  entropy  flux  relations. 
The  Onsagerist  approach,  on  the  other  hand,  is  basically  an  extension  of  classical 
thermodynamics  in  which  temperature  and  entropy  are  considered  as  physical  entities. 
On  this  basis,  second  law  considerations  lead  to  a  Clausius-Duhem  inequality  which  is 
different  from  that  postulated  by  the  rational  thermodynamicists. 

Recently,  criticisms  have  been  raised  against  the  rational  thermodynamic  approach 
by  Green  and  Naghdi  (1977,  1978),  Day  (1977)  and  Naghdi  (1980).  In  essence,  the 
criticisms  stem  from  the  inability  of  the  rational  thermodynamics  version  of  the 
Clausius-Duhem  inequality  to  rule  out  physically  unacceptable  behavior  in  some 
materials  and  an  unclear  relation  between  the  second  law  and  the  Clausius-Duhem 
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inequality.  Furthermore,  as  Day  (1977)  has  noted,  the  rational  thermodynamic 
approach  affords  no  way  of  giving  a  unique  meaning  to  the  concept  of  entropy.  As 
an  example,  consider  the  case  of  a  rigid  heat  conductor  in  thermal  equilibrium.  If 
heat  is  added  to  the  conductor,  the  form  of  the  Clausius-Duhem  inequality  that 
follows  from  rational  thermodynamics  predicts  that  the  subsequent  spatially 
homogeneous  temperature  field  of  the  conductor  decreases,  which  is  contradictory  to 
physical  experience  (Green  and  Naghdi,  1977).  The  Onsagerist  approach,  on  the  other 
hand,  does  not  suffer  from  the  above  criticisms  and  is  generally  accepted  by  the 
majority  of  researchers  in  the  field.  Accordingly,  we  adopt  the  Onsagerist  approach 
as  the  basis  for  the  present  study.  In  the  following  section,  the  pertinent 
thermodynamic  relations  that  follow  from  this  approach  are  presented. 

2.2  BASIC  RELATIONS  FROM  ONSAGERIST  APPROACH. 

It  is  well-known  that  the  state  of  a  reversible  system  can  be  completely 
described  by  functions  which  depend  only  on  observable  state  variables.  An 
irr.versible  system,  on  the  other  hand,  is  more  complicated  and,  to  define  its  state, 
requires  state  functions  that  depend  not  only  on  observable  state  variables  but  also  on 
internal  (non-observable)  variables  £  which  describe  the  internal  changes  that  occur  in 
a  material  during  an  irreversible  process.  Reversible  processes  can  then  be  treated 
within  the  framework  of  irreversible  thermodynamics  as  processes  for  which  the 
internal  variables  £  remain  fixed.  In  this  section,  the  basic  relations  of  the  Onsagerist 
school  of  irreversible  thermodynamics  with  internal  variables  are  presented;  the 
discussion  is  modeled  on  that  given  by  Nemat- Nasser  (1975.1978). 

The  first  law  of  thermodynamics  takes  the  following  rate  form 

5  =  e  -  2  •  e  (2.1) 

for  incremental  processes,  regardless  of  whether  they  are  reversible  or  irreversible.  Hn 
(2.1)  (  )  denotes  an  appropriate  objective  rate  of  ().  Here,  Q  is  the  heat  added  per 
unit  mass,  e  denotes  the  internal  energy  per  unit  mass,  £  represents  the  stress 
tensor,  and  £  is  the  small  strain  tensor. 

The  second  law  of  thermodynamics  requires  the  definition  of  entropy  and  an 
expression  for  entropy  production.  For  the  case  of  reversible  systems,  the  existence  of 
entropy  was  first  proven  by  Caratheodory  (see  Buchdahl,  1966).  Shortly  thereafter. 
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Valanis  (1971),  using  a  physical  argument,  addressed  the  question  of  entropy  for 
irreversible  systems  and,  in  a  more  recent  work  (1983),  he  formally  proved  the 
existence  of  entropy  for  such  systems  mathematically.  Using  the  concept  of  partial 
integrability  of  Pfaffian  forms,  Valanis  finds  the  temperature  6  =  0(g,e,£)  >  0  as  an 
integrating  factor  of  Eq.  (2.1),  which  renders  $/0  =  »?  as  an  exact  differential  form 
for  L  ~  L0  (=  constant),  where  r\  is  called  the  entropy: 

V  =  7U,e,j£)  (2.2) 


For  the  case  g  =  g  .  the  differential  form  of  Eq.  (2.2)  and  Eq.  (2.1)  yields 


•  £ 


i 


(2.3) 


from  which  it  follows  that 


2  = 


1-02 
0  ~  9e 


(2.4) 


For  an  arbitrary  irreversible  process  for  which  £  t  0,  Eq.  (2.2)  together  with 
Eqs.  (2.4)  leads  to  the  expression: 


*02 


e  ♦ 


§2 

9e 


:  4  2a 

e  *C 


(2.5) 


02 


t 


Upon  decomposing  ^  into  a  reversible  part  and  a  dissipative  part  Eq.  (2.5)  can 
be  rewritten  in  the  form: 


0tj  =  9rir  +  0*r}6 

where 


(2.6) 


(2.7) 


I 
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The  second  law  of  thermodynamics  states  that  the  rate  of  irreversible  entropy 
production  rj  must  be  either  zero  (for  reversible  processes)  or  positive  (for  irreversible 
processes).  Thus,  we  must  have 


(2.8) 


Alternately,  the  state  variables  e  and  rj  can  be  expressed  in  terms  of  £,  9  and  £ 
as  follows: 


e  =  .  V  =  =  V[i.6.£) 


(2.9) 


Then,  for  the  case  in  which  £  -  £q.  Eq.  (2.5)  can  be  rewritten  as 


*  .  •  +  92  1  _  1 

^  de  ~  *  90  ~  6 


#  •  0e  * 

ag  ~  ae  9 


i 

e 


9e 

H  "  1 


From  Eq.  (2.10),  it  follows  that 


(2.10) 


02  1 

a£  =  e 


9e  _ 
a£  '  2 


(2.11a) 


02  _  1  0e 
Q9  =  6  90 


(2.11b) 


The  Helmholtz  free  energy  f  is  now  introduced,  where 


£)  5  •  -  *7 

From  Eq.  (2.11a)  it  follows  that 


(2.12) 


Furthermore,  by  taking  the  derivative  of  Eq.  (2.12)  with  respect  to  6  and  using 
Eq.  (2.11b),  one  obtains  the  result 


V  = 


a* 

a  e 


(2.14) 


In  order  to  express  Eq.  (2.8)  in  terms  of  the  variables,  e.  6  and  £.  let  us  write 


i  (*,»,£)  =  | j  ■>(£.«<£.»,£).£) 


_  3a  «.  3a 
"a*  9; 


3e 

n 


Upon  combining  Eqws.  (2.4b)  and  (2.15),  it  folows  that 


(2.15) 


a  92  _  _  9e  +  M  _  _  9 1 
9  H  ’  9£  ■  9^ 


(2.16) 


As  a  result,  the  second  law  inequality  (2.8)  can  be  placed  in  the  form 

5  =  A  •  l  >  0  (2.17) 

where  the  internal  force,  is  defined  by 

4  -  -  ||  (218) 

In  a  formulation  based  upon  taking  £,  0  and  £  as  the  independent  variables, 
Eq.  (2.13)  may  be  solved  to  give 

£  •  £(£.«.£)  •  (2») 

Upon  introducing  the  Legendre  transformation: 

*  =  2  :  £  -  1  .  (2-2°) 
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it  follows  from  Eq.  (2.20)  that 


_  M 

~  “  9  a  ' 

M 

(2.21) 

9$ 

V  ~  90  ' 

(2.22) 

.  9$ 

k  *  u  ■ 

(2.23) 

The  rate-form  of  Eq.  (2.21).  with  the 

strain-rate  decomposition: 

1  =  ie  ♦  ip  ♦  lg  > 

(2.24) 

leads  to  the  following  relations 

•  _  } 

~e  9g9g  *•  ’  ~p  9g9£  * 

•  92* 

~  9g9 6 

(2.25) 

where  the  subscripts  e,  p.  and  0.  respectively,  refer  to  the  elastic,  plastic  and  thermal 
components  of  the  total  strain  rate. 

Note  that  in  Eq.  (2.25a),  92$/9g9g  is  a  tangent  compliance  matrix. 

2.3  CONDITIONS  FOR  THE  EXISTENCE  OF  A  PLASTIC  POTENTIAL. 

Non-associated  plasticity  rests  on  the  basic  premise  that  the  plastic  strain  rate 
is  defined  in  terms  of  a  scalar  plastic  potential  function  l  as  follows: 


where  the  scalar  X  is  known  as  the  plastic  multiplier.  We  will  now  explore,  within 
the  framework  of  the  irreversible  thermodynamics  just  described,  the  conditions  under 
which  such  a  relationship  exists. 

The  evolution  equations  for  the  internal  variables  £  are  commonly  taken  in  the 
general  form: 


l  «  a  (A.q 


(2.27) 


which  implies  that  the  rate  at  which  any  specific  structural  rearrangement  occurs  is 
fully  determined  by  the  current  values  of  the  associated  generalized  force  4  and  £. 
The  influence  of  the  macroscopic  stress  state  on  £  appears  only  implicitly  through  the 
dependence  of  A  on  a.  which  can  only  be  determined  from  micromechanical 
considerations  relating  the  generalizeJ  force  4  to  the  macroscopic  stress  a  (see  Rice. 
1971). 

In  order  to  prove  the  existence  of  a  plastic  potential  function  it  is  necessary 
to  assume  that  the  function  H  is  derivable  from  a  scalar  potential  0(4.£).  such  that: 


l  =  H  =  X 


M 


where 


A 


M) 


'4  B  •  dl 

X 


(2.28) 


(2.29) 


and  X  is  an  integrating  denominator  which  is  chosen  to  be  positive,  i.e., 

*(iU)  >0  for  . 

From  Eq.  (2.25),  we  can  write 

?  -  &L  l 

~p  9g9£  *  ' 


which  can  be  combined  with  Eq.  (2.28)  to  yield 


i  -  x  -  44  *  X 
*p  ‘  ae  84  - 


a 

92 


(2.30) 


(2.31) 


(2 . 32) 


since 
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(2 . 33) 


(2.34) 


* 

With  a  change  of  variable,  we  can  write 


i(A.£)  --  *[°-L) 


so  that  Eq.  (2.32)  can  be  written  as 


e 

~P 


=  X 


9^ 
9  a 


(2.35) 


which  is  identical  to  Eq.  (2.26)  and  shows  that  ^(g,£)  is  the  plastic  flow  potential 
with  X  the  plastic  multiplier.  To  obtain  this  expression,  it  was  necessary  to  assume 
that  Eq.  (2.28)  holds.  In  that  which  follows,  we  examine  the  conditions  under  which 
this  assumption  is  valid. 

It  has  been  pointed  out  (Nemat-Nasser,  1978)  that,  due  to  the  second  law, 
expressed  by  Eq.  (2.17),  \  cannot  be  prescribed  arbitrarily  in  the  same  manner  as  g 
or  £.  Also,  the  second  law  does  not  restrict  the  sign  of  the  term  Based  on 

the  sign  of  this  term,  a  stability  concept  has  been  introduced  (Valanis,  1968,  1971c: 
Nemat-Nasser,  1978).  A  material  is  said  to  be  thermodynamically  stable  if: 


l  •  l  <  0 


(2.36) 


and  otherwise  is  thermodynamically  unstable.  For  stable  materials.  j\  relaxes  when 
£  =  0,  while  for  unstable  materials,  ft  increases  when  \  —  0.  On  the  basis  of 
physical  considerations,  thermodynamically  unstable  materials  are  excluded  from  further 
discussion. 

Restricting  attention  now  to  thermodynamically  stable  materials,  let  us  consider 
the  conditions  that  must  be  imposed  on  the  function  H  to  insure  the  existence  of  a 
potential  0  (or  0).  Specifically,  it  can  be  shown  (see  Nemat-Nasser,  1975a:  1978) 
that  $  (or  0)  will  exist  if  either  one  of  the  following  two  conditions  are  satisfied: 

(i)  The  generalized  Onsager  relation  holds 


*  Again,  this  relationship  can  be  defined  explicitly  only  on  the  basis  of 
micromechanical  considerations. 


** 


This  concept  of  stability  is  not  related  to  the  stability  postulate  of  Drucker 
(1954). 
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(2.37) 


a 

fH . 

_L 

a 

fH.) 

i 

3A . 

lx 

■  3A . 

.X 

J 


where  X  is  an  integrating  denominator; 

(ii)  the  function  H  depends  only  on  A,  i.e.. 

H  =  H(&)  ,  (2.38) 


and  the  number  of  internal  variables  f.  is  less  than  8. 

*i 

Condition  (i)  is  satisfied,  for  example,  when  H  depends  linearly  on  X,  i.e., 
H.  =  B. .  A.  , 

i  ij  j  * 


(2.39) 


where  B..  represents  a  symmetric  positive  definite  matrix  to  satisfy  the  second  law 
inequality  (2.17).  Condition  (ii),  on  the  other  hand,  is  satisfied  when 


u  = 


(2.40) 


Again,  the  validity  of  the  above  assumptions  can  be  judged  only  on  the  basis  of 
micromechanical  considerations  of  the  response  of  the  particular  material  under 
consideration  which  is  outside  of  the  scope  of  this  report. 

2.4  RESTRICTIONS  IMPOSED  BY  THE  SECOND  LAW. 

There  are  basically  two  situations  for  which  strong  statements  can  be  made 
about  the  restrictions  imposed  by  the  second  law  on  non-associated  flow  rules.  In 
the  first  situation,  the  potential  function  ^  must  be  a  homogeneous  function  of  the 
generalized  force  4-  while  in  the  second  situation  the  Helmholtz  free  energy  function  f 
must  be  prescribed,  together  with  evolutionary  equations  for  the  internal  variables. 

Considering  now  the  first  situation,  we  can  write  from  Eqs.  (2.17)  and  (2.28) 
with  (2.34)  that 
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(2.41) 


5  =  4  .  X  ||  =  nX*(J,£)  =  nX*(2,£)  >  0 

N 


Since  n  >  0,  the  second  law  requires  only  that  X0S(g,£)  >  0.  If.  as  is  customary,  ^  is 
taken  to  be  positive,  the  the  second  law  requires  that 

X  >  0  ,  (2.42) 


which  is  consistent  with  Eq.  (2.30)  given  earlier.  Thus,  for  a  positive  potential 
function  l  which  is  homogeneous  in  fa,  Eq.  (2.42)  expresses  the  only  restriction 
imposed  by  the  second  law.  Inequality  (2.42)  is  often  used  to  distinguish  plastic 
loading  from  unloading  when  the  stress  state  is  on  the  yield  surface. 

Consider  now  a  special  subclass  of  the  above  materials  for  which  ^  is  a 
homogeneous  function  of  and  £  is  related  to  g  through  a  relation 

t  *  4(2  -  £„.£]  .  (2  43) 

such  that 

*[4(£  -  £„.£)]  -  r(e  -  «„,£).  (2.44) 

and  l  is  a  homogeneous  function  of  degree  m  (>o)  of  the  argument  g  -  g  .  Here, 
g  is  called  the  back  stress  and  generally  depends  on  £.  This  subclass  includes  a 
wide  variety  of  engineering  materials  whose  plastic  yield  characteristics  exhibit  both 
kinematic  and  isotropic  hardening.  For  these  materials,  the  plastic  potential  can  be 
expressed  in  the  form: 


r(4.£)  » r(e.t)  =  r0[e  -  £„(£)]  *  ^[£]  ■  (2  *s) 

Here,  f  describes  kinematic  hardening.  Since  ^  is  assumed  to  be  a 
homogeneous  function  of  degree  m(>0)  of  the  argument  g  -  go(£).  the  same  will  be 
true  for  ^  .  With  this  in  mind,  it  follows  that 
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e  =  X  =  X  f—  <p  fa  -  g  "I 
~p  og  og  7ol~  *-oJ 


(2.46) 


As  a  result,  the  rate  of  plastic  work.  W  .  can  be  extressed  as 

P 


\  =  (2  -  Zo) 


(2.47) 


so  that,  for  a  positive  plastic  potential  and  X  >  0,  it  follows  that 

ft  >  0  for  e  t  0 
P  ~P 


(2.48) 


Thus,  the  second  law  requires  positive  plastic  work  for  those  materials  whose  plastic 
potentials  are  of  the  form  of  Eq.  (2.45). 

In  the  second  situation  mentioned  above,  the  Helmholtz  free  energy  function  f 
and  the  evolutionary  equation  for  £  are  prescribed.  In  this  case,  the  second  law.  as 
reflected  by  Eq.  (2.17).  imposes  the  following  condition: 

-  ff  1  >-  0  ,  (2.49) 


which  places  restrictions  on  the  forms  that  f  and  £  may  take.  This  approach  was 
taken  by  Valanis  (1971)  in  developing  the  endochronic  theory.  For  this  purpose,  he 
postulated  the  simplest  isotropic  forms  for  the  Helmholtz  function  f  and  the 
evolutionary  equations  for  £  which  satisfy  Eq.  (2.49).  With  these,  he  was  able  to 
integrate  the  evolutionary  equations  and  subsequently  eliminate  the  internal  variables 
from  the  stress-strain  relations.  In  later  work,  Valanis  (1980)  showed  that  this 
approach  can  be  used  to  derive  the  equations  of  classical  associated  plasticity  with 
both  kinematic  and  isotropic  hardening. 

In  summary,  it  follows  from  the  discussion  above  that,  unless  the  Helmholtz  free 
energy  function  f  and  the  evolutionary  equations  for  £  are  known,  the  only  restriction 
imposed  by  the  second  law  on  non-associated  flow  rules  is  that  the  plastic  multiplier 
X  must  be  positive.  Since  the  general  approach  followed  in  practice  is  to  postulate 
forms  for  the  plastic  potential  function  f  and  the  yield  function  f,  and  not  to 
prescribe  either  f  or  £,  the  restriction  imposed  by  the  second  law,  X  >  0,  is  very 
weak  and  offers  little  constraint  to  the  non-associated  flow  rule. 
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SECTION  3 


THE  POSTULATES  OF  DRUCKER  AND  IL’IUSHIN 

In  the  past,  the  postulates  of  Drucker  (1960)  and  ll'iushin  (1961)  have  been 
adopted  as  bases  for  developing  certain  plasticity  theories  and  for  assessing  the 
acceptability  of  others,  despite  the  fact  that  these  postulates  are  not  based  on  the 
laws  of  thermodynamics.  The  purpose  of  this  section  is  to  show  that  these 
postulates  can  be  violated  without  causing  a  violation  of  the  second  law  of 
thermodynamics.  As  a  result,  it  follows  that  they  should  not  be  used  as  criteria  to 
judge  the  thermodynamic  acceptability  of  inelastic  constitutive  theories. 

3.1  DRUCKER’S  POSTULATES. 

Consider  an  elasto-plastic  material  in  equilibrium  under  a  state  of  stress  g°  and 
£°.  denoted  by  state  0.  An  external  agency  slowly  applies  an  independent  set  of 
forces  in  equilibrium,  reaching  state  1  with  stress  g'  and  strain  £ ' .  Drucker’s 
postulate  of  stability  in  the  large  requires  that  the  following  inequality  hold: 

V 

[o  -  g°)  •  de  >  0  (3.1) 

o 

h 

Over  a  cycle  of  application  and  removal  of  the  additional  stress,  it  is  also  postulated 
that  the  net  work  performed  by  the  external  agency  is  non-negative: 

t  * 

e 

~P 

[o  -  2°]  ‘  dep  >  0  ,  (3.2) 


where  £°  and  e'  denote  the  plastic  strains  at  states  0  and  1.  respectively.  If  these 
states  are  an  infinitesima1  distance  apart,  we  can  write  g’  -  g°  =  Ag  and  Eqs.  (3.1) 
and  (3.2)  reduce  to  the  expressions: 
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La  •  Ae  >  0  , 

*4  e*  ' 


(3.3) 


La  •  Le  >0 

~  p 


From  Eq.  (2.24)  in  the  preceding  section,  we  have  for  isothermal  processes: 

.  -  jll  A/t  AC 

~p  3g9cr  ^  +  d£3g  ^ 


Ae  =  Ae  ♦  Ae  = 

fv  — 


(3.4) 


Since  the  elastic  compliance  is  positive  definite,  it  follows  that 


g2* 

*2  •  a2e  *  4«  •  3^ 


fig  >  0  for  fig  *  g 


(3.5) 


Inasmuch  as  the  inequalities  (3.3b)  and  (3.5)  are  sufficient  conditions  for  Eq.  (3.3a), 
only  the  inequality  (3.3b)  needs  to  be  considered.  Using  Eq.  (3.4)  we  write 


&2  •  %  =  a2  •  44  ^ 


(3.6) 


which,  in  view  of  Eqs.  (2.17)  and  (2.23),  can  be  rewritten  in  the  form: 
La  •  Ae  =  Ag  •  ^  C 


(3.7) 


This  equation  expresses  Drucker's  postulate  of  stability  in  the  small  in  terms  of  the 
increment  of  dissipation,  AD.  Since  the  second  Law  requires  only  that  AD  >  0,  it 
follows  that  the  inequality  (3.7)  can  be  violated  for  processes  in  which  the  second  law 
is  not  violated. 

3.2  ILMUSHIN'S  POSTULATE. 

A  postulate  introduced  by  ll’iushin  (1961)  states  that  if  the  work  AW  done  by 
external  forces  in  an  isothermal  closed  cycle  of  strain  is  positive,  plastic  deformation 
takes  place  and,  if  the  work  AW  is  zero,  only  elastic  deformation  occurs.  Expressed 
mathematically,  this  postulate  takes  the  form: 
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AW  =  6  a  •  de  >  0  , 

T  M  M  ’ 


(3.8) 


where  the  integration  is  taken  around  a  closed  path  in  strain  space.  Since  a  strain 
cycle  is  always  possible,  even  for  softenin^stress-strain  behavior,  ll’iushin’s  postulate 
is  less  restrictive  than  Drucker’s  postulate. 


Consider  a  small  strain  cycle,  such  that  eQ  *  +  A£  £2  =  £  .  during 

which  the  stress  changes  as  follows:  £  ■*£  •♦  a The  work  done  during  tfiis  cycle 

is  given  by  the  expression: 


(VA£ 


AW  = 


2  *  <*£  ♦ 


2  *  d£  =  2  K  *  2l) 


Ae  -  f  (2l 


22)  *  Ae  ■  (3.9) 


e  +Ae 


Hence 

tw  =  5  (20  -  «2)  •  *2 


(3  10) 


Let  us  now  assume  that  irreversible  changes  occur  only  during  loading.  Then  the 
internal  variables  £  will  evolve  as  follows: 

f0  -  £,  =  £„  ♦  '  £i  ■  (3.ii) 

From  Eq.  (2.13),  we  can  write 

20  -  22  =  af  '  4  ^°,£°  *  '  (312) 


Drucker’s  postulate  expressed  by  the  inequality  (3.2)  requires  a  stress  cycle. 
Such  a  cycle  is  not  possible  on  the  softening  branch  for  materials  with 
softening  stress-strain  relations. 
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so  that 


o 

~0 


22 


-  .  -O 


(3.13) 


Thus,  in  view  of  Eqs.  (3.8),  (3.10).  and  (3.13),  we  can  write 

af  (»)  •  Hi  °  , 


(3.14) 


which  expresses  ll’iushin’s  postulate  in  terms  of  the  gradient  of  the  dissipation 
increment.  9(AD)/9e.  and  the  amplitude  of  the  strain  cycle,  Ae.  Since  the  second  law 
requires  only  that  AD  >  0,  it  follows  that  the  inequality  (3.14)  can  be  violated 
without  violating  the  second  law.  Thus,  the  work  done  in  a  closed  strain  cycle  may 
be  negative  without  violating  the  laws  of  thermodynamics. 

3.3  CLOSURE. 

In  this  section,  we  have  expressed  both  Drucker’s  and  ll’iushin’s  postulates  in 
terms  of  the  dissipation  increment,  AD.  so  that  a  direct  connection  with  the  second 
law  can  be  explored.  This  appears  to  be  a  new  result.  An  examination  of  the 
resulting  expressions  clearly  shows  that  both  postulates  are  not  related  to  the  second 
law  of  thermodynamics. 
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SECTION  4 


N ON-ASSOCIATED  FLOW  LAWS 

In  the  preceding  section,  the  restrictions  imposed  by  thermodynamics  on  elasto- 
plastic  constitutive  theories,  in  general,  were  established;  these  restrictions  apply  to 
both  associated  and  non-associated  flow  laws.  In  this  section,  attention  is  specifically 
focused  on  non-associated  flow  laws,  and  the  implications  of  the  thermodynamic 
restrictions  for  such  laws  are  explored.  In  addition,  an  unusual  feature  of  non- 
associated  plasticity  is  discussed,  and  it  is  shown  that  non-associated  flow  laws 
satisfy  the  continuity  condition  of  Prager  (1949).  Finally,  path-dependence  of  non- 
associated  plasticity  is  considered. 

4.1  RESTRICTION  IMPOSED  BY  THE  SECOND  LAW. 

Consider  an  elastic-plastic  material  having  a  yield  criterion  of  the  form: 

f (e.*)  =  o  (4.1) 


where  x  is  some  scalar  hardening  parameter.  The  elastic  portion  of  the  deformation 
is  governed  by  an  incremental  form  of  Hooke’s  Law,  which  we  write  as: 


e  =  -  £P3 


(4-2) 


where  L  is  the  elastic  modulus  tensor.  The  plastic  strain  rate  £  is  taken  to  follow 
the  non-associated  flow  law: 


=  X 


U 

92 


(4.3) 


where  ^  denotes  the  plastic  potential  function.  The  difficulty  of  not  being  able  to 
eliminate  £  in  Eq.  (2.30)  in  terms  of  observable  variables  for  \  t  Q  due  to  the  lack 
of  the  information  on  the  microstructural  irreversible  processes  has  been  circumvented 
by  the  use  of  Prager’s  consistency  condition.  The  consistency  condition  for  the  yield 
surface  f(g,x),  with  plastic  work  hardening  given  by  x  =  g  •  leads  to  the  result 


18 


ii  . 
92 


1=  •  i 

A# 

M  _  If  TM 

3g  ”  3x  *  gg 


(4.4) 


It  is  straightforward  to  show  that  the  numerator  of  Eq.  (4.4)  is  positive  for  plastic 
loading  and  negative  for  elastic  unloading.  Therefore,  in  order  to  satisfy  the  general 
restriction  imposed  by  the  second  law  of  thermodynamics,  namely,  X  >  0.  it  is 
necessary  that 


3f 

92 


9f 

3g  ”  3  *2 


(4.5) 


since  the  denominator  of  Eq.  (4.4)  is  a  function  of  g  and  *  and  not  of  g  and  g. 
The  inequality  (4.5)  provides  the  basic  restriction  imposed  on  non-associated  plastic 
flow  laws  by  thermodynamics  for  materials  which  exhibit  plastic  work  hardening. 
When  there  is  no  hardening,  i.e„  perfect  plasticity,  the  inequality  (4.5)  reduces  to  the 
condition 


3f 

92 


So 

A# 


(4.6) 


which  was  first  proposed  by  Bleich  (1972).  We  shall  return  to  Eq.  (4.6)  later  on  in 
Section  6  to  explore  specific  restrictions  on  constitutive  models  of  particular  interest. 

4.2  THE  CONTINUITY  CONDITION. 

The  continuity  condition,  introduced  by  Prager  (1949),  is  a  general  requirement 
for  elasto-plastic  constitutive  theories  which  is  designed  to  prohibit  possible  difficulties 
that  might  arise  from  the  presence  of  awkward  discontinuities  in  constitutive  relations 
for  neutral  or  near-neutral  loading.  Difficulties  in  using  constitutive  relations  which 
violate  the  continuity  condition  have  been  reported  in  the  past  (Nelson  and  Baladi, 
1977).  In  essence,  the  continuity  condition  requires  that  a  neutral  loading  condition 
be  a  limiting  case  of  either  loading  or  unloading.  As  a  result,  the  constitutive 
relations  for  loading  and  unloading  must  coincide  for  this  limiting  case  so  that  there 
are  no  discontinuities  in  the  constitutive  response.  The  purpose  of  this  section  is  to 
demonstrate  that  non-associated  plastic  flow  laws  do,  in  fact,  satisfy  the  continuity 
condition. 
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Consider  a  work  hardening  material  having  a  yield  function  in  the  form  of 
Eq.  (4.1),  with 


K  - 


2  •  de 

N  IV  v 


(4.7) 


From  the  consistency  condition  and  Eq.  (4.3),  it  follows  that 

)  *°  • 


IH*X 


E(«-E 


(4  6) 


from  which 


X  = 


0f  • 

* 


5*  54  5g 


(4.9) 


subject  to  the  requirement  that 

-  SJ  «  •  if  >  0  ■ 


(4.10) 


Let  us  now  rewrite  Eq.  4.9)  as 

^H(e.«)||2||(!!f  -n„) 


(4.11) 


where 
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h(2  *)  = 


.  Si  ,  .  M 

dii  °  do 


(4.12) 


Substitution  of  Eq.  (4.11)  into  Eq.  (4.3)  leads  to  the  result 


£p  =  H(2'*}l|2|l(Cf  *  fla)  Be 


(4.13) 


where  is  defined  by  replacing  f  of  (4.12b)  with 
Taking  the  norm  of  Eq.  (4.13),  we  obtain 


=  cose 


(4.14) 


where 


flf  *  %  =  cos0 


(4.15) 


Equation  (4.14)  defines  the  manner  in  which  the  magnitude  of  the  plastic  strain  rate 
vector,  ll|  II,  varies  with  the  direction  and  magnitude  of  the  stress  rate  vector. 
From  this  expression,  we  see  that 


I  im  • 

e+*rf 2  ~p  ~ 


(4.16) 


Hence,  the  non-associated  flow  law  satisfies  the  continuity  condition. 
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4.3  A  PECULIAR  FEATURE  OF  NON-ASSOCIATED  FLOW  LAWS. 


Non-associated  flow  laws  exhibit  an  unusual  feature  of  response  that  is  not 
shared  by  associated  flow  laws.  In  particular,  for  a  range  of  loading  directions  in 
stress  space,  the  non-associated  flow  law  can  lead  to  response  which  is  stiffer  than 
purely  elastic  response.  In  this  section,  we  examine  this  feature  of  non-associated 
plasticity,  following  the  approach  originally  taken  by  Mroz  (1966). 

Consider  an  arbitrary  stress  rate  vector,  g,  and  let  e  denote  the  projection  of 
,  •  •  .  **  o 

the  strain  rate  vector,  £.  on  g.  i.e.. 


(4.17) 


Define  directional  and  elastic  moduli  as  follows: 


e 


a 


so  that  from  Eqs.  (4.17)  and  (4.18),  we  can  write 


K 


e 


Z 


e 

~e 


(4.18) 


(4.19) 


From  the  strain  rate  decomposition 


e  = 


e 

~e 


♦  e 
~P 


< 


(4.20) 


it  follows  that 


•  •  •  •  •  • 

£  •  i  =  Z  *  £e  ♦  Z  *  £p  • 

Let  us  rewrite  Eq.  (4.13)  in  the  following  form 
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(4.21) 


a 


cosd 


(4.22) 


1 1  i 


I 


where 


K  5 


(4.23) 


and 


c°s 8  =  J3f  *  na 

Substitution  of  Eqs.  (4.19)  and  (4.22)  into  Eq.  (4.21)  results  is 


From  Figure  1,  we  can  write 
=  cos(0  ♦  p) 


(4.24) 


(4.25) 


Thus.  Eq.  (4.24)  can  be  placed  in  the  form: 

K-  =  K~  +  K  cos(^  +  P)cos&  (4.26) 

a  e 

Note  that  when  0  <  151  <  t/2  and  ir/2  <  15  +  ^1  <  r,  it  follows  from  Eq.  (4.26) 
that 


Kd>Ke  .  (4.27) 

A  typical  non-dimensional  plot  of  Eq.  (4.26)  is  given  in  Figure  2.  where  the 
dimensionless  directional  modulus.  K  ./K  ,  is  shown  as  a  function  of  6  for  different 

O'  € 

values  of  /?.  As  this  figure  reveals,  the  modulus  K  .  can  exceed  the  corresponding 
elastic  modulus  for  a  range  of  values  of  p.  This  feature  of  non-associated 
plasticity  is  not  shared  by  associated  flow  laws. 
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Kd/Ke 


Figure  1  Geometry  of  the  potential  and  yield  surfaces  at  a  point  in 
stress  space,  showing  the  angles  6  and  p. 


Figure  2. 


Typical  features  of  Eq.  (4.26)  for  various  values  of  the 
parameter  p. 


4.4  PATH  DEPENDENCE. 

Another  response  feature  of  non-associated  plasticity  which  is  not  found  in 
associated  plasticity  is  path-dependence.  To  show  path  dependence  of  non-associated 
plasticity,  consider  the  following  system  of  equations  which  define  a  general  non- 
associated  elasto-plastic  material: 


2  =  1=  fi  -  ip]  . 

(4.28) 

*P  _ 

(4.29) 

N 

II 

H- 

(4.30) 

Here,  Eq.  (4.28)  is  the  incremental  form  of  Hooke’s  Law,  with  [.  the  elastic  modulus 
tensor:  Eq.  (4.29)  is  the  non-associated  flow  law,  with  l  being  the  plastic  potential. 
In  addition,  f  denotes  the  yield  function,  which  depends  on  g  and  a  scalar  hardening 
parameter,  k,  where  *  =  g  •  |p. 

From  the  consistency  condition  for  f(g,x)  =  0,  it  follows  that 


9f 

a2 


2 


(4.31) 


By  combining  the  above  equations,  it  can  be  shown  that 


0f 


a2 


•  l  •  i 

M 

WToTTTu 

9g  “  9*  2  g g 


(4.32) 


which  was  given  earlier  in  Eq.  (4.4).  Upon  introducing  Eq.  (4.32)  into  Eq.  (4.29), 
and  substituting  the  result  into  Eq.  (4.28),  the  following  expression  is  obtained: 


A  comparison  of  Eqs.  (4.38)  and  (4.39)  shows  that,  when  t  f,  we  have 


BkCij  *  Bi  j  k£. 


(4.40) 


Consider  now  a  rate-potential  $/  defined  such  that 


If  such  a  potential  exists,  it  follows  that 


02ft 


a2fl 


3a 


k£. 


9a.  . 
i  J 


3e,0  0e. 0de..  0e..deLO  3e.. 

k<L  ki.  i  j  i  j  kIL  i  j 


which,  in  view  of  Eq.  (4.34),  implies  that 


BijkC  =  Bk£ij  ' 


(4.41) 


(4.42) 


(4.43) 


Thus,  since  Eq.  (4.43)  holds  for  associated  flow  rules,  a  rate  potential  exists  in  this 
case.  On  the  other  hand,  in  view  of  Eq.  (4.40),  Eq.  (4.42)  is  violated  by  non- 
associated  flow  rules  and  therefore  a  rate-potential  does  not  exist  for  this  case.  This 
knowledge  will  be  useful  in  the  developments  which  follow. 


Let  us  now  consider  the  path-dependency  of  the  strain  energy  function.  For  this 

purpose,  we  consider  an  incremental  process  from  time  t  to  t  +  At,  and  denote  the 

0,0 

stress  and  strain  at  time  t  by  a  and  e  .  Over  the  time  interval  At,  the  increments 
in  stress  and  strain  are  then  given  by  the  expressions: 


A2l 


’t  +At 

o 

Jdt 

t 

o 


rfco*At 


(4.44) 


The  change  in  the  strain  energy  during  the  incremental  process  is: 


27 


AW  = 


o  •  edt  = 

N  A* 


edt 


(4.45) 


ft^At 


t 

o 


ft^+At 


fWAt 


edt  ♦  a 

~  ~c 


t 

O 


which  can  be  placed  in  the  form 


fAe. 
~1 


AW  = 


A  c 


*  <*[H]  ♦  20  *  Hi  . 


(4.46) 


where 

rt 

edt  for  tQ  <  t  <  t  ♦  At  .  (4.47) 

t 

o 


Ag(t)  =  2(t)  -  eQ  ,  Ae  = 


Consider  two  different  paths  in  Ae-space.  as  depicted  in  Figure  3.  namely.  I  :  Q 

*  H  *  H-  and  ^2:  ^  *  H  *  Hi-  F°r  eac^  of  these  paths,  let  us  evaluate  the 
first  integral  on  the  right-hand  side  of  Eq.  (4.46): 


(4.48) 


aw2  = 


A o  •  dfAe 


(ii) 


Then 


AWj  -  AW2  = 


2) 


H  •  <*{H)  =  J  ii_i2a2  •  d(A0 


(4.49) 


(4.50) 


where  the  path  1^  -  1^  implies  Q  ■»  A£  ■*  A^^  ■»  H  +  Q.  Upon  applying  Stoke's 
Theorem  in  the  non-dimensional  strain  increment  space  to  transform  the  line  integral 
to  a  surface  integral  (see  Flanders,  1963),  one  obtains: 
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AW 


1 


h-h 


>  A o  ■  .  d  Cac  '  -0 

(b>  ,j 


z: 

S  ( 1 » j )  <  (k,«0 


f9A<r 


k£ 


9A  a  ■  -1 
_ LL 


9Ae 


'J 


9Ae 


k£ 


6{Le.i}]K6[^ 


(4.51) 


when  S  is  chosen  so  that  the  boundary  of  S:9S  =  I  -  l2  and  A  denotes  the  exterior 
product,  i.e.. 


d (ac  1  -0  Ad Cac k^0  =  -  dfoe^Adfoe'-O  (4.52) 

and  where  (i.j)  <  (k.£.)  implies  that  the  pair  (k,£.)  follows  the  pair  (i.j)  in  an 
established  sequence  of  pairs. 

In  view  of  Eq.  (4.40).  the  integrand  of  the  integral  in  Eq.  (4.51)  is  nonzero  for 
non-associated  flow  rules  and  zero  for  associated  flow  rules.  Therefore,  for  non- 
associated  flow  rules  we  may  have 

AWj  -  AW2  t  0  (4.53) 


while,  for  associated  flow  rules, 
AW 1  -  AW2  =  0 


(4.54) 


This  proves  the  path-dependency  of  the  strain  energy  integral.  It  is  important  to  note 
that  the  strain  energy  integral  exists  uniquely  even  for  non-associated  flow  rules  since 
there  always  exists  one  and  only  one  g  for  a  given  e,  as  shown  in  Eq.  (4.34). 

In  the  incremental  numerical  approaches  used  in  finite  element  and  finite 

difference  methods,  it  is  tacitly  assumed  that  Ae  ,  is  achieved  in  a  proportional 

manner,  that  is,  for  any  t  such  that  t  <  t  <  t  +  At,  we  have 

o  o 

Ae(t)  =  Afij  . 
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(4.55) 


The  implication  of  the  path-dependency  of  AW  for  non-associated  flow  rules,  therefore, 
is  that  in  the  numerical  approaches  one  must  use  small  time  steps  so  that  the 
proportional  path  approximation  does  not  significantly  deviate  from  the  actual  strain 
path.  For  associated  flow  rules,  there  is  no  concern  since  the  strain  path  does  not 
change  the  integral  AW  for  a  specified  terminal  strain  increment,  Ae. 

The  above  requirement  for  small  time  steps  in  the  integration  of  the  equations 
for  non-associated  plasticity  is  physical  and  independent  of  other  considerations  based 
upon  numerical  stability. 
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SECTION  5 


UNIQUENESS 

In  this  section,  the  question  of  uniqueness  of  initial-boundary  value  problems  for 
elasto-plastic  materials  having  non-associated  flow  rules  is  explored.  Uniqueness  is,  of 
course,  an  issue  that  is  completely  unrelated  to  thermodynamics.  Valanis  (1985) 
showed  that  the  usual  approach  for  proving  uniqueness  for  associated  elasto-plastic 
materials  is  valid  under  the  assumption  of  proportional  paths  even  when  the  plastic 
flow  is  non-associated.  However,  when  the  requirement  is  not  satisfied,  the 
bifurcation  and  stability  analyses  require  a  special  care  for  non-associated  plasticity 
theories. 

5.1  SUFFICIENT  CONDITION. 

Consider  a  three-dimensional  continuum  at  time  t  having  volume  V  and  surface 
S.  Let  g(x,t),  £(x.t)  and  v(x.t)  denote  the  stress  field,  strain  field  and  velocity  field, 
respectively,  in  V.  In  addition,  let  g(x,0)  =  Q,  £(x.O)  =  Q  and  Ji(x.O)  =  0  in  V. 
Also,  let  ST  be  that  part  of  the  surface  over  which  the  traction  vector  X(t)  is 
prescribed  and  SD  the  remaining  portion  over  which  the  displacement  is  prescribed  in 
the  time  domain  0  <  t  <  «. 

Restricting  attention  to  small  deformations,  we  proceed  to  explore  the  conditions 
under  which  the  solution  g(xt),  £(x,t),  and  ^(x,t)  in  V  is  unique,  for  0  <  t  <  *, 
when  the  following  equations  apply: 

In  V: 


a 


'  j  *  j 


(5.1) 


a  ■  . 

'  J 


=  BijkCekC 


(5.2) 


where  (5.1)  is  the  equation  of  motion  in  the  absence  of  a  body  force  and  (5.2)  is  the 
basic  constitutive  equation  for  an  elasto-plastic  solid.  For  a  non-associated  flow  rule, 
the  inelastic  tangent  modulus,  B..k^,  is  non-symmetric,  i.e.. 


Bi  jkfc.  *  Bk£.i j 


(5.3) 
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a  fact  that  will  prove  to  be  important  in  the  developments  which  follow. 

On  ST: 

a.  .n.  =  T.  ;  (5.4) 

ij  j  i 


On  S 


D' 


v.  is  prescribed 


(5.5) 


Consider  now  an  incremental  process  from  t^  to  =  t  +  At.  and  assume  that 

a  unique  solution  exists  up  to  time  tQ.  In  order  to  prove  uniqueness  by  negation,  we 

assume  two  solutions  1  and  2  of  the  incremental  process  and  define  the  difference 

state  between  t  and  t„: 

o  1 


Av . 


I 


.CD 


A  a  ■ 


=  oW  -  aW 


J  • 


J  ' 


J  ' 


(5.6) 


for  which  the  values  at  and  tj  are  prescribed: 

i  (*■«)  *  ‘"j ;  Co)  * 0  ■ 

**i2)W  -  *il)W  •  (5'7) 

JOj  .  (tj  i  ^  i  .  ojf  (t,)  -  (tj  . 

Since  the  two  solution  states  satisfy  (5.1),  the  difference  state  satisfies 

A<7 ■  •  -  flAv.  =  0  (5-8) 

J',J  r  ' 


Multiplying  Eq.  (5.8)  by  (-Av.)  and  integrating  over  the  domain  V  and  time  between 

t  and  t,,  one  obtains 
o  1 
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I  = 


-ho..  .Id  V  dt 
J '  ,JJ 


o 


lit 


Pv;)2}  - 


♦  ho ■  ■ Av  .  . 
J  >  '  »J 


dVdt  =  0 


(5.9) 


Using  the  Gauss  theorem  and  the  boundary  conditions  for  the  difference  state  obtained 
from  (5.4)  and  (5.5).  Eq.  (5.9)  reduces  to 


(5.10) 


(5.11) 


The  first  integral  of  (5.10)  defines  the  kinetic  energy  and  is  zero  only  when  Av.  =  0. 
For  the  associated  flow  problems  with  an  incremental  potential  AA/ 


SAW 


SAe .  . 
•J 


(5.12) 


since  the  stress  increment  ho.,  is  a  unique  function  of  he...  Introducing  AW  into 
(5.10),  one  obtains 


34 


'  r 

I  =  {|  (a7  .  ]  2  ♦  AwjdV  =  0 


(5.13) 


in  which  the  kinetic  energy  of  the  difference  state  and  A/V  are  both  positive  definite. 
Therefore,  in  order  to  satisfy  (5.10),  we  must  have: 


Av.  =0  and  Ac..  =0  everywhere  in  V. 

i  i  j 


This  is  contrary  to  the  assumption  that  Av.  =  0  and  Ac..  ^  0  in  V.  Therefore, 
uniqueness  has  been  proven  by  negation. 

The  above  procedure  works  well  for  associated  flow  plasticity  in  which  a 
potential  exists  and  Eq.  (5.12)  holds. 

The  approach  requires  special  care,  however,  for  non-associated  flow  rules.  Since 
in  this  case  A  a.,  is  path-dependent,  it  is  possible  to  have 

AV.  =  Ac..  =0  and  yet  to.  ■  ^  0  (5. 14) 

i  ij  ’  ij 

for  special  path. 

However,  restricting  to  proportional  paths  in  achieving  Ae..  and  to..,  the 
inequality 

AW  =  J  to. .  Ae. .  >  0  (5.15) 

2  i  j  i  j 

furnishes  a  sufficient  condition  for  uniqueness  but  not  a  necessary  condition  even 
under  quasi-static  load  environment  (Valanis,  1985). 

It  should  be  noted  that  for  sufficiently  small  At,  the  o..  -  e..  path  followed  to 
achieve  Ae..  in  an  incremental  process  may  be  approximated  by  a  piecewise 
incremental  proportional  path.  Adopting  this  assumption,  Mroz  (1963:1966)  developed 
a  special  class  of  non-associated  flow  models  based  upon  the  uniqueness  requirement 
so  that  they  possess  the  minimum  hardening  required  to  satisfy  the  conditions  (5.15). 
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If  Ac..  is  decomposed  into  its  elastic  and  plastic  components,  the  condition  (5.15) 
can  then  be  written  in  the  form 


Aa.  .Ae. .  =  A«7. .  [Ac.  .]  ♦  Atf.  •  [Ae.  .1  >0  . 

U  ij  ijl  'Jt  U'P 


(5.16) 


The  first  itim  on  the  right-hand  side  of  (5.16)  denotes  the  incremental  elastic  strain 
energy  and  is  positive  definite.  The  second  term,  associated  with  the  plastic 
deformation,  becomes  negative  only  in  the  case  of  non-associated  flow  laws  for  stress 
increments  A  a.,  whose  direction  6  falls  within  the  range  t/2  <  \9  +  ()\  <  *  (see 
Section  4.3).  J 

Hill  (1958)  has  shown  that  the  requirements  (5.15)  is  a  sufficient,  but  not 
necessary,  condition  for  stability.  Also,  the  violation  of  this  condition  is  known  to 
cause  some  types  of  localized  failure.  The  judgement  of  whether  or  not  ^to  accept 
such  local  failures  as  real  or  not  should  be  based  on  experimental  evidence. 

5.2  DISCUSSION. 

In  this  section,  the  question  of  uniqueness  for  initial  value  problems  involving 
non-associated  elasto-plastic  materials  has  been  explored.  We  have  shown  that  the 
usual  mathematical  approach  for  proving  uniqueness  of  initial  value  problems  involving 
elasto-plastic  solids  with  associated  flow  laws  holds  even  for  solids  with  non- 
associated  flow  laws.  In  several  cases,  plasticity  theories  of  non-associated  flow  laws 
have  been  constructed  based  upon  the  uniqueness  requirement,  such  as  for  example 
those  considered  by  Mroz  (1963,  1966).  Parenthetically,  we  note  that  the  inability  to 
prove  uniqueness  should  not  be  taken  to  imply  that  uniqueness  does  not  exist:  it 
simply  means  that  we  are  unable  to  answer  the  question  of  uniqueness.  To  address 
uniqueness  of  non-associated  flow  problems,  a  new  mathematical  approach  is  needed, 
which  is  beyond  the  scope  of  this  study.  At  this  moment,  it  is  interesting  to  note 
the  inadequacy  of  associated  flow  plasticity  in  predicting  localized  failure,  such  as 
necking,  in  metals.  For  example,  the  numerical’^tudy  of  necking  in  plane  strain 
tensile  test  by  McMeeking  and  Rice  (1975)  which  employed  a  associated  flow 
theory  found  continued  growth  of  the  diffuse  neck,  with  no  tendency  toward  intense 
localized  shearing.  Successful  simulations  by  Tvergaard,  Needleman  and  Lo  (1981), 
and  Iwakuma  and  Nemat-Nasser  (1982)  utilized  the  flow  theory  which  includes  non- 


For  further  discussion  of  destabilizing  effects  of  non-associated  flow  laws,  see 
Maier  and  Hueckel  (1979). 
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coaxiality.  For  pressure-sensitive  geomaterials.  Rudnicki  and  Rice  (1975)  utilized  the 
nonuniqueness  of  non-associated  flow  theories  to  predict  the  localization  of  deformation 
as  a  constitutive  instability  (bifurcation).  The  above  examples  indicate  that  the  non- 
associated  flow  rules  which  exhibit  non-normality  and  non-coaxiaility  of  stress  tensor 
and  plastic  strain  rate  tensor  may  be  of  necessity  to  predict  the  inception  of  unstable 
deformation  at  realistic  stress  level. 
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SECTION  6 


APPLICATION  TO  SOME  NON-ASSOCIATED  MODELS 
HAVING  .INFLOW  POTENTIALS 

In  Section  4.1,  the  restrictions  imposed  on  non-associated  flow  models  by 
thermodynamics  were  explored.  It  was  shown  that  the  second  law  of  thermodynamics 
is  satisfied  whenever  the  inequality  given  earlier  in  Eq.  (4.6).  i.e.. 


9f 
5  a 


(6.1) 


is  satisfied.  The  issue  of  uniqueness,  discussed  in  Section  5,  was  much  less 
conclusive. 

In  this  section,  we  consider  two  non-associated  flow  models  and  investigate  the 
restrictions,  if  any.  imposed  upon  them  by  the  second  law  requirement  (6.1).  The 
specific  models  were  selected  for  examination  because  they  represent  two  important 
classes  of  models  presently  being  used  by  the  defense  community  in  ground  motion 
simulation  studies.  Both  models  are  elastic-perfectly  plastic,  and  thus  do  not  exhibit 
strain  hardening.  Also,  they  both  have  pressure-dependent  yield  surfaces  with 
J2-plastic  potentials. 

6.1  GENERALIZED  DRUCKER-PRAGER  MODEL. 

This  model  has  a  pressure  dependent  yield  surface  of  the  form: 
f  (*iP)  =  J2  -  F(p)  ,  (6.2) 


and  a  von  Mises-type  plastic  potential,  i.e. 

m  =  j2 


(6.3) 


where  g  denotes  the  deviatoric  stress  tensor,  p  is  the  pressure,  and 


•j 


(6.4) 
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The  function  F  is  positive  and  increases  monotonically  with  p. 

The  deviatoric  response  follows  an  incremental  form  of  Hooke’s  Law 


where  e  and  ep  denote,  respectively,  the  deviatoric  strain  tensor  and  the  plastic 
deviatoric  strain  tensor,  while  fi  denotes  the  shear  modulus.  The  hydrostatic  behavior 
follows  a  hypoelastic  expression  of  the  form: 


P  =  K(e)e  (6.6) 

where  p  is  the  pressure,  e  represents  the  volumetric  strain,  and  K  denotes  the 
tangent  bulk  modulus. 

From  Eqs.  (6.5)  and  (6.6).  L  in  Eq.  (4.2)  may  be  expressed  as: 

UijM>  *  H  -  I  V  •  <6'7> 

Upon  substituting  Eqs.  (6.2),  (6.3)  and  (6.7)  into  the  left-hand  side  of  Eq.  (6.1).  one 
obtains  the  result: 


8f  |  9 

3ff .  .  i  jk li  do.  0 
ij  k  t 


4^2 


(6.8) 


Thus,  since.  J2  >  0,  the  model  satisfies  the  condition  (6.1)  without  any  restrictions. 

6.2  A  THIRD  INVARIANT  GENERALIZED  DRUCKER-PRAGER  MODEL. 

This  model  is  similar  to  that  described  in  the  preceding  section  except  that  the 
yield  surface  is  of  the  form: 

f(§.P]  =  J2  '  fjp)  J3  '  p2(P)  =  0  (6-9) 


where  F(p)  is  a  positive,  monotonically  increasing  function  of  p,  0  <  a  <  1,  and 
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(6.10) 

The  flow  potential  is  given  by  Eq.  (6.3).  and  Eqs.  (6.5)  and  (6.6)  also  apply,  with  L 
defined  by  Eq.  (6.7).  It  is  then  straightforward  to  show  that  the  left-hand  side  of 
Eq.  (6.1)  takes  the  form: 


I  A 

••  0S  ■  •  S  .  |  s, 

3  3  ij  jk  ki 


3f 

9a .  . 
U 


Lijk£ 


2J, 


(6.11) 


To  prove  the  Eq.  (6.11)  is  positive  for  0  <  a  <  1,  let  us  introduce  the 
variables  r ,9  in  the  f-plane  as  follows: 


r 


2 


♦ 


cos30  = 


3J3 

2 


(6.12) 

(6.13) 


Equation  (6.9)  can  then  be  rewritten  as: 


f  _  d.  . 

2  F(p) 


— —  cos3 6  -  F^(p)  =  0 


(6.14) 


On  f  =  0  in  the  ir-plane,  r  is  strictly  a  function  of  9  and 

jk 

at  sin  36  =  0 

Also,  for  cos30  =  0,  it  follows  from  Eq.  (6.14)  that 


as  depicted  in  Figure  4. 


(6.15) 


(6.16) 
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Equation  (6.11)  may  be  rewritten,  using  Eqs.  (6.12)  and  (6.13)  as: 


For  cos30  <  0,  it  is  trivial  to  show  that  I  >  0.  To  prove  that  I  >  0  for 
cos30  >  0,  which  includes  0  =  0,  +  2r/3  at  which  r  takes  maximum  values  (see 
Eq.  (6.15)),  first  note  that 


(f(p)]  *  for  0  -  a  <  1 


(6.18) 


max 


which,  when  used  in  Eq.  (6.18),  leads  to  the  result: 


I  >  2(i2{\  -  — 

1  ft 


F(p)Jmax 


}  *  2Mr2(l -a)  >  0 


(6.19) 


Thus,  the  model  satisfies  the  basic  thermodynamic  requirement,  Eq.  (6.1).  without  any 
restrictions. 
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SECTION  7 


CONCLUSIONS 


In  this  study,  we  have  explored  the  possible  restrictions  that  thermodynamics 
and  uniqueness  impose  on  problems  involving  non-associated  elasto-plastic  solids. 
From  this  investigation,  the  following  major  conclusions  have  emerged: 


1.  In  the  case  of  elasto-plastic  constitutive  models  formulated  by 
defining  a  plastic  potential  a  yield  surface  f  and  a  non-associated 
flow  rule,  i.e., 


the  only  restriction  imposed  by  thermodynamics  is  that  the  plastic 
multiplier  \  must  be  positive,  i.e.,  X  >  0. 

2.  If  an  inelastic  constitutive  model  is  constructed  by  defining  the 
Helmholtz  free  energy  function  f  and  the  evolutionary  equations  for 
the  internal  variables,  £,  thermodynamic  considerations  then  require 
that  the  following  condition  be  satisfied: 


& 


i  *  o 


This  condition  usually  leads  to  restrictions  on  the  form  of  the 
evolutionary  equations  for  the  internal  variables 

3.  The  postulates  of  Drucker  and  ll'iushin  are  shown  to  be  unrelated 
to  the  second  law  of  thermodynamics.  It  was  found  that  these 
postulates  can  be  violated  for  processes  in  which  the  second  law  is 
satisfied. 

4.  Non-associated  flow  laws  satisfy  the  continuity  condition  of  Prager 
(1949). 
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5. 


Elasto-plastic  solids  having  non-associated  flow  laws  exhibit  a 
seemingly  peculiar  response  feature  in  which  the  directional  plastic 
modulus  can  exceed  the  corresponding  elastic  modulus  for  a  certain 
range  of  stress  increment  vectors  in  stress  space.  For  this 
response,  nonunique  solutions  may  exist,  which  suggest  the 
inception  of  localized  failure  due  to  bifurcation. 

6.  One  of  the  major  differences  between  associated  anc  non-associated 
flow  laws  is  that  associated  flow  laws  are  not  path  dependent  while 
non-associated  flow  laws  are  path-dependent. 

7.  Because  of  the  path-dependent  nature  of  non-associated  flow  laws, 
the  approach  for  proving  uniqueness  for  initial  value  problems 
involving  associated  flow  laws  remains  valid  for  exploring  uniqueness 
under  the  proportional-path  assumption  even  if  the  flow  law  is  non- 
associated. 

8.  Uniqueness  has  been  proven  only  for  a  special  class  of  non- 
associated  flow  laws  in  which  there  is  sufficient  hardening  to 
eliminate  destabilizing  effects.  For  non-associated  model  with  no 
strain  hardening,  uniqueness  cannot  be  proven  by  following  the 
usual  approach  without  performing  bifurcation  and  stability  analyses. 
A  new  mathematical  approach  is  needed  for  exploring  uniqueness  of 
non-associated  flow  problems. 

9.  Two  generic  types  of  non-associated  elasto-plastic  models  that  are 
commonly  used  by  the  defense  community  to  simulate  explosive 
ground  motion  were  considered  in  detail.  It  was  shown  that  both 
types  of  models  satisfy  the  second  law  as  is,  without  the  need  for 
further  restriction.  On  the  other  hand,  since  both  models  are 
perfectly  plastic  (no  hardening),  there  is  a  potential  for  instability 
and  possibly  a  lack  of  uniqueness.  Also,  care  must  be  taken  in 
numerically  dealing  with  non-associated  flow  rules  since  they  lead  to 
non-selfadjoint  operators,  which  are  known  to  require  special 
treatment. 

I 
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On  the  basis  of  the  above  remarks,  it  follows  that  the  greatest  concerns  in 
using  non-associated  flow  models  without  hardening  arise  from  (1)  the  uncertainty  of 
the  question  of  uniqueness,  (b)  possible  instabilities,  (c)  numerical  difficulties  in 
dealing  with  non-selfadjoint  operators  and  (d)  the  seemingly  peculiar  feature  of 
response  in  which  the  directional  plastic  modulus  can  exceed  the  corresponding  elastic 
modulus  when  instability  is  anticipated.  Further  research  is  clearly  needed  to  resolve 
these  issues  before  non-associated  flow  laws  can  be  used  with  reasonable  confidence. 
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JOINT  DATA  SYSTEM  SUPPORT  CTR 
ATTN:  R  MASON 
ATTN:  C-332 

JOINT  STRAT  TGT  PLANNING  STAFF 
ATTN:  JK 
ATTN:  JKCS 
ATTN:  JLT 
ATTN:  JP 
ATTN:  JPEP 

LAWRENCE  LIVERMORE  NATIONAL  LABORATORY 
ATTN:  DNA-LL 

NATIONAL  DEFENSE  UNIVERSITY 
ATTN:  ICAF-ICC 
ATTN:  NWCO 

NATIONAL  SECURITY  AGENCY 
ATTN:  CHIEF  A  GROUP 

OFFICE  OF  THE  JOINT  CHIEFS  OF  STAFF 
ATTN:  ED30 
ATTN:  J-3 
ATTN:  J-5 


3  CYS  ATTN:  J-5  NUCLEAR  DIV/STRATEGY  DIV 
ATTN:  J-8/CAD 
ATTN:  JAD/SFD 
ATTN.  JAD/SSD 

OFFICE  OF  THE  SEC  OF  DEFENSE 
ATTN.  NAVAL  FORCES 
2  CYS  ATTN:  STRAT  PROGRAMS  &  TNF 

STRATEGIC  AND  THEATER  NUCLEAR  FORCES 
ATTN:  DR  SCHNEITER 

U  S  NATIONAL  MILITARY  REPRESENTATIVE 
ATTN:  U  S  DOCUMENTS  OFFICER 

UNDER  SEC  OF  DEFENSE  FOR  POLICY 
ATTN.  DUSP/P 
ATTN:  USD/P 

UNDER  SECRETARY  OF  DEFENSE 
ATTN:  DEPUNDSEC 
ATTN:  GSEVIN 

DEPARTMENT  OF  THE  ARMY 

DEP  CH  OF  STAFF  FOR  OPS  &  PLANS 
ATTN:  DAMO-NCN 
ATTN:  DAMO  NCN 
ATTN:  DAMO-RQS 
ATTN:  DAMO-ZXA 

DEP  CH  OF  STAFF  FOR  RSCH  DEV  &  ACQ 
ATTN:  DAMA-CSM-N 

HARRY  DIAMOND  LABORATORIES 
ATTN:  SLCIS-IM  TL 

U  S  ARMY  BALLISTIC  RESEARCH  LAB 
ATTN:  SLCBR  DD  T 
ATTN:  SLCBR-SS-T 

U  S  ARMY  COMB  ARMS  COMBAT  DEV  ACTY 
ATTN:  ATZL-CAP 

U  S  ARMY  COMD  &  GENERAL  STAFF  COLLEGE 
ATTN.  ACQ  LIBRARY  DIV 
ATTN:  ATZL  SWT  A 

U  S  ARMY  MATERIEL  COMMAND 
ATTN.  AMCCN 

U  S  ARMY  MATERIEL  SYS  ANALYSIS  ACTVY 
ATTN:  AMXSY-CR 

U  S  ARMY  MISSILE  COMMAND 
ATTN:  AMSMI-XF 

U  S  ARMY  NUCLEAR  &  CHEMICAL  AGENCY 
ATTN:  MONA-NU 

U  S  ARMY  NUCLEAR  EFFECTS  LABORATORY 
ATTN:  DR  J  MEASON 

U  S  ARMY  WAR  COLLEGE 
ATTN:  LIBRARY 
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DEPARTMENT  OF  THE  NAVY 

NAVAL  INTELLIGENCE  SUPPORT  CTR 
ATTN  NISC-30 

NAVAL  OCEAN  SYSTEMS  CENTER 
ATTN:  CODE  9642 

NAVAL  POSTGRADUATE  SCHOOL 
ATTN:  COOL  1424 

NAVAL  RESEARCH  LABORATORY 
ATTN:  CODE  1220 
ATTN:  CODE  2627 

NAVAL  SEA  SYSTEMS  COMMAND 
ATTN:  SEA  09G53 
ATTN:  SEA  643 

NAVAL  WEAPONS  EVALUATION  FACILITY 
ATTN:  CLASSIFIED  LIBRARY 

NUCLEAR  WEAPONS  TNG  GROUP,  ATLANTIC 
ATTN:  CODE  221 
ATTN:  DOCUMENT  CONTROL 

NUCLEAR  WEAPONS  TNG  GROUP,  PACIFIC 
ATTN:  DOCUMENT  CONTROL 

OFC  OF  THE  DEPUTY  CHIEF  OF  NAVAL  OPS 
ATTN:  OP  654 
ATTN:  OP  981 

OFFICE  OF  THE  CHIEF  OF  NAVAL  OPERATIONS 
ATTN:  CNO  EXECUTIVE  PANEL 

OPERATIONAL  TEST  &  EVALUATION  FORCE 
ATTN:  CODE  80 

OPERATIONAL  TEST  &  EVALUATION  FORCE, 
ATTN:  INTEL  OFFICER 

TACTICAL  TRAINING  GROUP,  PAC'FIC 
ATTN:  COMMANDER 

THEATER  NUCLEAR  WARFARE  PROGRAM  OFC 
ATTN:  PMS423 

DEPARTMENT  OF  THE  AIR  FORCE 

AFIA/INKD 

ATTN:  MAJ  COOK 

AFIS/INT 

ATTN:  INT 

AIR  FORCE  SYSTEMS  COMMAND 
ATTN:  DL 
ATTN:  SD 

AIR  FORCE  WEAPONS  LABORATORY 
ATTN:  SUL 

AIR  UNIVERSITY  LIBRARY 
ATTN:  AUL  LSE 


ASSISTANT  CHIEF  OF  THE  AIR  FORCE 
ATTN:  SAF/ALR 

DEPUTY  CHIEF  OF  STAFF/AF  RDQM 
ATTN:  AF/RDQI 

SPACE  DIVISIGN/CWH 
ATTN:  CWH 

STRATEGIC  AIR  COMMAND/NRI  STINFO 
ATTN:  NRI/STINFO 

STRATEGIC  AIR  COMMAND/STIC 
ATTN:  STIC 

DEPARTMENT  OF  ENERGY 

LAWRENCE  LIVERMORE  NATIONAL  LAB 
ATTN:  L-35 
ATTN:  L-38 
ATTN:  L-389 
ATTN:  W  HOGAN 
ATTN:  L-53 

ATTN:  Z  DIVISION  LIBRARY 

LOS  ALAMOS  NATIONAL  LABORATORY 
ATTN:  D  RICHMOND 
ATTN:  TDOWLER 
ATTN:  REPORT  LIBRARY 

SANDIA  NATIONAL  LABORATORIES 
ATTN:  TECH  LIB  3141 
ATTN:  RB  STRATTON 

OTHER  GOVERNMENT 

CENTRAL  INTELLIGENCE  AGENCY 

ATTN:  N10-  STRATEGIC  SYS 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

ADVANCED  RESEARCH  &  APPLICATIONS  CORP 
ATTN:  DOCUMENT  CONTROL 

BDM  CORP 

ATTN:  CWASAFF 
ATTN:  J  BODE 
ATTN:  J  BRADDOCK 
ATTN:  R BUCHANAN 

BOEING  CO 

ATTN:  H  WICKLEIN 
ATTN:  JW  RUSSELL 

COMPUTER  SCIENCES  CORPORATION 
ATTN:  F  EISENBARTH 

GRUMMAN  DATA  SYSTEMS  CORPORATION 
ATTN:  SSHRIER 

HORIZONS  TECHNOLOGY,  INC 
ATTN:  J  PALMER 

I  IT  RESEARCH  INSTITUTE 

ATTN:  DOCUMENTS  LIBRARY 
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INSTITUTE  FOR  DEFENSE  ANALYSES 
ATTN:  CLASSIFIED  LIBRARY 
ATTN:  JGROTE 

JAYCOR 

ATTN:  R  SULLIVAN 

KAMAN  SCIENCES  CORP 
ATTN:  F  SHELTON 

KAMAN  SCIENCES  CORP 
ATTN:  E CONRAD 

KAMAN  SCIENCES  CORPORATION 
ATTN.  DASIAC 

KAMAN  SCIENCES  CORPORATION 
ATTN:  DASIAC 

MARTIN  MARIETTA  DENVER  AEROSPACE 
ATTN  J  DONATHAN 


ORION  RESEARCH  INC 

ATTN:  JESCHOLZ 

PACIFIC-SIERRA  RESEARCH  CORP 
ATTN  H  BRODE 

PACIFIC-SIERRA  RESEARCH  CORP 
ATTN:  D  GORMLEY 


R  &  D  ASSOCIATES 

ATTN:  C  K  B  LEE 
ATTN:  0  SIMONS 

2  CYS  ATTN:  DOCUMENT  CONTROL 

R  &  D  ASSOCIATES 

ATTN:  C  KNOWLES 
ATTN:  J  THOMPSON 


R  &  0  ASSOCIATES 

ATTN.  GGANONG 


RAND  CORP 

ATTN:  V  JACKSON 


RAND  CORP 

ATTN.  B BENNETT 


S-CUBED 

ATTN 
2 CYS  ATTN 
2 CYS  ATTN 


B  PYATT 
H  MURAKAMI 
H  READ 


SCIENCE  APPLICATIONS  INTL  CORP 
ATTN.  DOCUMENT  CONTROL 
ATTN:  ESWICK 
ATTN:  J  MARTIN 
ATTN:  M  DRAKE 
ATTN:  R  J  BEYSTER 


SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  J  SHANNON 
ATTN.  LGOURE 
ATTN:  WLAYSON 


SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  RCRAVER 

TRW  INC 

ATTN:  D  SC  ALLY 
ATTN:  R  BURNETT 


> 


Dist-3 


